Absence of correlations in the QCD Dirac spectrum at high temperature 
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I propose a simple model of the distribution of the small eigenvalues of the QCD Dirac operator 
well above the finite temperature phase transition where chiral symmetry is restored and the spectral 
density at zero vanishes. Assuming the absence of correlations between different regions of the low 
lying spectrum I derive analytic formulas for the distribution of the first two eigenvalues. I find 
good agreement with data obtained using the overlap Dirac operator in quenched SU{2) lattice 
simulations. This suggests that if chiral symmetry is restored spectral correlations are not important 
and all the statistical properties of the spectrum are encoded in the spectral density. 
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Random Matrix Theory (RMT) has been rather suc- 
cessful in describing some universal properties of the 
spectrum of the QCD Dirac operator Most of the 
results, however, concern the chirally broken phase of the 
theory. This regime is characterized by a non-zero den- 
sity of eigenvalues around zero, which through the Banks- 
Casher relation leads to a non-vanishing chiral conden- 
sate. Using a single parameter, the value of the chiral 
condensate, RMT makes detailed predictions about the 
distribution of the smallest Dirac eigenvalues. 

Above the finite temperature phase transition, T^, the 
Dirac spectrum is much less understood in terms of RMT. 
Multicritical random matrices with fine tuning in the ac- 
tion might be able to describe the chirally restored phase 
0. Another possibility is to add a temperature depen- 
dent constant matrix to the usual chiral random matrix 
to mimic the lowest Matsubara mode 0, 0| ■ For compar- 
isons with lattice staggered data above Tc see Refs. Q. 
Recently there has also been renewed interest in a better 
understanding of the Dirac spectrum around and above 
Tc with the hope of clarifying the connection between the 
chiral and deconfinement transition 

In the present paper I propose a model of the low lying 
Dirac spectrum that can provide a simple alternative to 
RMT well above Tc . The main assumptions I make for 
the Dirac spectrum around zero are the following: 

(1) For fixed temperature the spectral density scales 
with the spatial volume. 

(2) The spectral density per unit spatial volume is 



p{\) - CA", 
where C and a are constants. 



(1) 



(3) The number of eigenvalues in any two disjoint in- 
tervals are independent random variables. 



TABLE I: The number of configurations A'^ generated for the 
different spatial box sizes L. Nq-q is the number of configu- 
rations in the trivial topological sector, the one used here. 



L 


12 


14 


16 


18 


20 


22 


24 


32 


N 


900 


850 


738 


400 


490 


326 


376 


42 


Nq=o 


879 


821 


711 


379 


451 


298 


328 


29 



Using lattice data from quenched SU{2) simulations with 
the overlap Dirac operator I directly verify assumptions 
(1) and (2). Based on (1-3) I derive analytic formulas for 
the statistical properties of the first two eigenvalues in 
terms of C, a and the spatial volume. Fitting C and a 
to the lattice data, the analytic formulas provide several 
parameter-free predictions and I find perfect agreement 
with the numerical lattice data. 

At first I summarize the details of the numerical sim- 
ulations. The data is based on quenched simulations of 
the SU{2) gauge theory with Wilson plaquette couphng 
f3 = 2.6 and time extension Nt = 4. This corresponds 
to a temperature of about T = 2.6Tc, well above the fi- 
nite temperature phase transition. The spatial size of the 
box was chosen in the range TV^ = 12 — 32, spanning more 
than an order of magnitude in spatial volume (see Table 
m. In eight different spatial volumes I computed the 16 
(or 32) eigenvalues of smallest magnitude of the overlap 
Dirac operator 8], 



Do^ = l^A [A^A] " 



A^l + S-Dn 



(2) 
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where Dq is the Wilson Dirac operator and I use s — 0.4 
that appears to be optimal for the condition number of 
A^ A. The spectrum of the Dirac operator is symmetric 
with respect to the real axis and I only consider the eigen- 
values with non-negative imaginary parts. The overlap 
Dirac operator has an exact chiral symmetry and its spec- 
trum lies along a circle touching the origin. The eigen- 
values of smallest magnitude therefore have only small 
real parts and the low end of the spectrum is much like 
that of the continuum Dirac operator which has a purely 
imaginary spectrum. In the analysis I always use only 
the imaginary part of the eigenvalues. 
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(a) (b) 




FIG. 1: The cumulative spectral density scaled by the 3- 
volume for all eight 3- volumes. The dashed line is the an- 
alytic prediction using the parameters fitted for the average 
smallest eigenvalue. The difference between Fig. (a) and (b) 
is that while in the former we assume no gap for the fitting, 
in the latter a gap G — 0.39 is assumed. 

The low-end of the Dirac spectrum is known to depend 
strongly on the temporal fermionic boundary condition 
which is effectively a combination of the Polyakov loop 
and the explicitly chosen anti-periodic boundary condi- 
tion 9] . In the quenched SU (2) theory the Polyakov loop 
Z{2) symmetry is spontaneously broken above the crit- 
ical temperature. While in the negative Polyakov loop 
sector there is a non-zero density of eigenvalues around 
zero (as if chiral symmetry were broken), in the posi- 
tive sector there seems to be a gap roughly controlled by 
the lowest Matsubara frequency Although in the 

quenched theory the two sectors are equivalent, I will re- 
strict my study to the latter, the one that survives in 
the high temperature phase in the presence of dynamical 
fermions since they explicitly break the Z{2) symmetry. 
The boundary conditions I use for the fermions are al- 
ways anti-periodic in the time direction and periodic in 
all the spatial directions. 

In addition to the small non-zero modes, due to its 
exact chiral symmetry, the overlap operator can also have 
exact zero modes, regardless of the boundary condition 
or the Polyakov loop sector. I use only configurations 
belonging to the trivial topological sector which is the 
most abundant in these simulations (see Table |T]). 

(1) Scaling of the spectral density — In Fig. [1] I plot 
the cumulative spectral density normalized by the spatial 
volume, 

r(A)= / dxp{x), (3) 
Jo 

for different spatial box sizes. r(A) is the average num- 
ber of eigenvalues per unit volume smaller than A. The 



densities in the different volumes agree perfectly except 
for the very low end where the data is sparse and statis- 
tical errors are large. The scaling of the spectral density 
with the 3-volume is a rather non-trivial property of the 
interacting theory. For free fermions the level spacing 
is inversely proportional to the linear box size since the 
Dirac equation is first order. That would imply a spectral 
density scaling with the linear size, not the volume. 

(2) Form of the spectral density — As can be seen in the 
log-log plot of r, the integrated spectral density versus A 
(Fig. [1^) , the data is well described by the simple power 
law of Eq. ID) therefore I will assume that in what follows. 

I emphasize that for the main result of the paper the 
particular form of the spectral density is immaterial. The 
whole forthcoming analysis can be carried out with any 
simple analytic form of the spectral density that describes 
the data accurately. In particular there are also expecta- 
tions based partly on RMT [s^ and lattice data 711 that at 
high temperature, well above Tc the spectrum develops 
a gap G around zero. In that case the spectral density 
would be expected to start off from G as 

p{X)=C{X-Gr. (4) 

Later I will also discuss this possibility. Potentially the 
accumulation of small non-zero eigenvalues seen in ^10| 
might also complicate the picture, but the temperature 
used here is high enough that these modes are completely 
absent, at least with the present statistics. 

(3) No spectral correlations — i.e. the number of eigen- 
values in any two disjoint intervals are independent ran- 
dom variables. This is the most non-trivial of the three 
assumptions and the demonstration of its validity is the 
main result of the paper. I will show that the resulting de- 
scription fits the numerical data rather accurately. The 
absence of correlations between different regions of the 
spectrum immensely simplifies any model since it means 
that all probabilistic information about the spectrum is 
encoded in the spectral density. All the properties of 
the low eigenvalues follow from a simple few-parameter 
description of the spectral density like Eq. ([T]) or Q. 

It is important to note that the absence of correlations 
between different regions of the spectrum is not related 
in any simple way to the correlations between individual 
eigenvalues that is often considered in random matrix 
models [T2j. In particular it does not imply the absence 
of correlations between individual eigenvalues. 

Distribution of the smallest eigenvalues — I now com- 
pute the distribution of the smallest two eigenvalues from 
assumptions (1-3). From the definition of the spectral 
density Vp{X)dX is the average number of eigenvalues in 
spatial volume V in an interval of length dX centered 
around A. If at fixed volume the length of the interval 
dA — > then the probability of having one eigenvalue in 
the interval is Vp{X)dX and the probability of no eigen- 
value is 1 — Vp{X)dX. This is because the probability of 
having more than one eigenvalue becomes negligible. 

Let us first compute the probability Pno{Xi., X2) that 
there is no eigenvalue in the interval [Ai, A2]. Since corre- 
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lations between the subintervals are ignored the probabil- 
ity of having no eigenvalue in an interval can be written 
as the product of probabilities of having no eigenvalue in 
any subinterval of its decomposition into small subinter- 
vals. 



Pno{M, A2) 



lim (1 - Vp{xi)Ax) (1 



Vp{x2)Ax) ... 
(1 - Vp{xn)Ax),{5) 



where xi = Ai, xn = A2 and Xk+i—Xk — Ax. Expanding 
the product the terms can be organized in powers of Ax 
and in the limit Ax — > the order n term goes to the 
following n-fold integral; 



Sn 



/ dxip{xi) / dx2p{x2) 



Ai 



dXnPiXn)- 



(6) 



Substituting the spectral density in Eq. ([T]) for p the in- 
tegrations can be explicitly carried out resulting in 



S„ 



i-CV) 



a+l\ 



n\{a + ly 

Finally the summation over n is easily done giving 

CV 



Pnoi^l, A2) 



exp 



1 



(7) 



(8) 



for the probability of no eigenvalue in the interval [Ai , A2] . 

Using again assumption (3) the probability of having 
the smallest eigenvalue around A is the product of two 
probabilities: having no eigenvalue in [0, A] and having 
one eigenvalue in [A, A-|-dA]. Thus the probability density 
of the smallest eigenvalue is 



pi(A) = exp 



CV 

a+1' 



A° 



CVX". 



(9) 



The average smallest eigenvalue is easily calculated from 
this as 



(Ai)=/ piix)xdx = {CVp)-^Til + p), (10) 

JQ 

where I introduced the notation p = [1 + a)^^. 

The distribution and average of the second smallest 
eigenvalue can be calculated in a similar fashion. Again 
the probability of having the second eigenvalue around 
A can be decomposed as a product: having the smallest 
eigenvalue at x, having no eigenvalues between x and A 
and having an eigenvalue around A. Since x (the occur- 
rence of the first eigenvalue) can be anywhere between 
and A, x has to be integrated and the probability density 
of the second eigenvalue is 



P2(A) 



dx pi{x) Pno{x, A) p(\) 



exp I -A 



a+l 



a 



1 



(11) 
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FIG. 2: The average smallest and second smallest eigenvalue 
as a function of the spatial volume in lattice units. The solid 
curve is a two-parameter fit of the form Eq. pop . The dashed 
curve is the fitted parameter-free prediction for the average 
second smallest eigenvalue. 



Finally we obtain 

{\2)^{CVp)-^'T{2 + p), 



(12) 



for the average second eigenvalue. The procedure can be 
easily continued for higher eigenvalues. 

Numerical results and tests — In the remainder I com- 
pare the analytic formulas and the data. There are many 
different ways to proceed; one has to fit the two parame- 
ters C and a in the spectral density using any of the ana- 
lytic formulas above and then the rest of the formulas are 
predictions that can be tested against the data. I choose 
to use Eq. (|10p , the scaling of the average smallest eigen- 
value with the 3- volume V for the fit. Using a fit range of 
y = 143 - 32^ produces a good fit with x^/d.o.f.= 0.89. 
The parameters obtained are C — 1.07(82) x 10^ and 
/i = ^ = 0.0327(12). In Fig. [2] I plotted the average 
smallest eigenvalue as a function of the 3-volume along 
with the fit. Also in this Figure I present the average 
second eigenvalue and the analytic formula of Eq. (fT2|) 
using the above parameters. Notice that this is now a 
prediction of the model with no further parameters to be 
fitted and it agrees quite well with the data. 

A more detailed prediction is the distribution of the 
smallest and second smallest eigenvalue for different spa- 
tial volumes. As an illustration, Fig. [3] shows this in a 
spatial volume of = 16'^ and again there is good agree- 
ment with the data. The picture is similar for the other 
volumes down to ^ = 14"^ below which even the distribu- 
tion of the smallest eigenvalue starts to be concentrated 
above the upper limit of the validity of the simple power- 
law of Eq. ID) for the spectral density. 

Is there a spectral gap ? — Coming back to the question 
of whether there is a spectral gap, the derivation of the 
distributions can be easily generalized for the spectral 
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FIG. 3: The distribution of the smaUest and second smaUest 
eigenvalues in a spatial volume of = 16^. The solid lines 
indicate the analytic prediction of the model (Eqs. ([9} and 

density of Eq. @ yielding 

(Ai> - [CViiyvil + ii) + G (13) 

for the average smallest eigenvalue. Unfortunately a 
three-parameter fit — including G — to the average small- 
est eigenvalue is not very meaningful since in a wide range 
of G almost equally good fits can be found for C and a. 
We can, however, assume that G is somewhat smaller 
than the smallest observed eigenvalue. In our simula- 
tion data for the eight different spatial volumes compris- 



ing altogether more than 4000 configurations the smallest 
eigenvalue was 0.3965. If a gap of G = 0.39 is assumed, 
the fit yields C = 5.4(3.6) x 10^ and = 0.1275(52) with 
X^/d.o.f.= 1.03. Comparing the predicted cumulative 
spectral density with the data on a log-log plot (Figs. 
[T^ and Wp) shows that the cumulative spectral density 
definitely favors G = over G = 0.39, however, a much 
smaller but still non-zero value of the gap cannot be ruled 
out based on the present data. 

Discussion — We proposed a description of the low- 
lying Dirac spectrum at high temperature based on the 
assumption that there are no correlations between differ- 
ent regions in the spectrum. Our model provides detailed 
analytic predictions for the distribution of the low-lying 
Dirac eigenvalues and perfect agreement is found with 
numerical lattice data. The absence of correlations is a 
peculiar property of the chirally symmetric phase and 
is not expected to hold if chiral symmetry is broken. 
It would be interesting on the one hand to get a more 
fundamental understanding of this property and on the 
other hand to see whether it is shared by a wider class 
of chirally symmetric systems. A good testing ground 
could be the Schwinger model with two flavors of mass- 
less fermions [l^ . The model proposed here might even- 
tually provide a full description of the Dirac spectrum 
in the chirally restored phase, just like Random Matrix 
Theory in the chirally broken phase. Together with RMT 
it might also open the possibility to a better understand- 
ing of the nature of the chiral phase transition of QCD 
and even the conformal phase of QCD-like theories with 
more fermionic degrees of freedom [isj . 
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